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Abstract.  Observations  of  quasi-invariant  m  -3  power  spectra  of  various  atmospheric 
perturbations  at  large  vertical  wavenumbers  m  have  received  considerable  theoretical 
attention.  Yet  some  other  observations,  most  notably  in.  the  stratosphere,  have 
revealed  significant  departures  from  this  shape  and  much  reduced  spectral  densities  at 
large  m .  Here  it  is  argued  that  these  changes  arise  when  a  spectrum  of  gravity  waves 
encounter  mean  wind  changes  which  cause  intrinsic  horizontal  phase  speeds  to 
increase.  Concomitant  changes  in  vertical  wavenumbers,  amplitudes,  and  shear 
variance  can  produce  changes  in  spectral  character  qualitatively  similar  to 
observations.  The  two  important  dynamical  processes  governing  the  spectral  response 
are  refraction  of  vertical  wavenumbers  and  conservation  of  Eliassen-Palm  flux  for 
nondissipating  waves,  processes  already  encapsulated  in  both  single-wave  and  spectral 
parameterizations  of  gravity-wave  processes.  Retrospective  application  of  these  ideas 
to  various  observations  in  the  literature  leads  to  successful  “prediction”  (based  on  the 
background  wind  profile)  of  either  quasi-invariant  m~ 3  spectra  or  attenuated  spectra  at 
large  m.  However,  while  explaining  their  occurrence,  the  precise  changes  in  spectral  shape 
predicted  by  simple  models  of  wind-shifted  spectra  often  differ  from  those  observed. 


Introduction 

Spectral  analysis  of  vertical  profiles  of  horizontal-velocity 
(«'),  relative-temperature  and  relative-density  fluctuations 
in  the  atmosphere  reveal  wavenumber  power  spectra  of  the 
form  aA(N)m~q  at  large  wavenumbers  m,  where  N  is  the 
background  Brunt-Vaisala  frequency  N,  A(N)  is  a  polariza¬ 
tion  factor  which  depends  on  the  atmospheric  variable  being 
measured  (e.g.,  A(N)  =  N2  for  u'  data),  and  a  is  a  constant. 
Early  data  suggested  q  ~  2.4  [Van  Zandt,  1982],  but  a  range 
of  later  measurements  suggested  near-constant  values  of  q  = 
3  and  a  ~  1/10-1/6  [e.g.,  Smith  et  al.,  1987;  Fritts  et  al.,  1988; 
Tsuda  et  al.,  1989;  Wu  and  Widdel,  1989,  1991], 

Since  then,  a  number  of  models  have  been  developed 
which  seek  to  explain  these  quasi-invariant  spectral  shapes. 
Currently,  they  can  be  grouped  into  three  main  categories 
[e.g.,  Hines,  1993a]:  (model  a)  saturation  models,  in  which 
the  spectral  shape  is  due  to  saturated  wave  amplitudes  [e.g., 
Dewan  and  Good,  1986;  Smith  et  al.,  1987];  (model  b) 
Doppler-spread  models,  in  which  the  shape  arises  from 
stochastic  wavenumber  spreading  of  large-m  waves  by  the 
velocity  oscillations  of  smaller-zn  waves  [Hines,  1991b, 
1993b];  and  (model  c)  diffusive  models,  in  which  nonlinear 
interactions  among  the  waves  act  diffusively  upon  the  wave 
ensemble  to  yield  an  m~ 3  spectral  shape  [e.g.,  Weinstock, 
1990;  Gardner,  1994;  Zhu,  1994].  These  classifications  are 
somewhat  arbitrary,  as  there  are  similarities  and  differences 
among  these  and  other  proposed  models  [e.g.,  Muller  et  al., 
1986;  Dewan,  1994;  Sato  and  Yamada,  1994], 

However,  recent  higher-resolution  measurements,  per¬ 
formed  throughout  the  year  at  different  locations  and  alti- 
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tudes,  have  on  occasion  revealed  greater  variations  in  the 
large-m  portion  of  the  spectrum  than  most  models  currently 
allow.  Most  notably,  the  spectra  at  large  m  in  the  strato¬ 
sphere  are  often  (although  not  always)  less  intense  (a  — 
1/60-1/20)  and  shallower  or  steeper  than  expected  (q  ^  3) 
[e.g.,  Shibata  et  al.,  1988;  Wilson  et  al.,  1990,  1991b;  Tsuda 
etal.,  1991, 1994b ;  Senft  et  al.,  1993;  Hostetler  and  Gardner, 
1994;  Mitchell  et  al.,  1994;  Whiteway  and  Carswell,  1994; 
Allen  and  Vincent,  1995].  These  features,  when  they  arise, 
tend  to  be  persistent  and  climatological  rather  than  transi¬ 
tory  [e.g.,  Tsuda  et  al.,  1991;  Wilson  et  al.,  1991b].  Some 
mesospheric  data  also  exhibit  spectral  shape  and  intensity 
variations  at  large  m  which  exceed  model  constraints,  al¬ 
though  here  the  form  and  occurrence  of  the  spectral  changes 
seem  more  random  and  less  climatologically  reproducible 
[e.g.,  Wu  and  Widdel,  1990;  Senft  and  Gardner,  1991; 
Collins  et  al.,  1994]. 

Some  theoretical  spectral  models  have  already  been  ap¬ 
plied  or  adapted  to  account  for  these  observations.  Using  the 
saturation  hypothesis,  some  have  argued  that  attenuated 
shallow  spectra  such  as  these  indicate  an  unsaturated  wave 
field  [e.g.,  Tsuda  et  al.,  1991],  while  others  have  argued  that 
they  may  reflect  dynamical  saturation  of  inertia  gravity 
waves  [e.g.,  Shibata  et  al.,  1988;  Wilson  et  al.,  1990].  Hines 
[1993a]  showed  that  similar-looking  spectral  variations  can 
arise  within  a  Doppler-spreading  framework  by  considering 
input  source  spectra  of  different  forms.  Similar  effects  are 
accommodated  within  the  wave-diffusive  model  of  Wein¬ 
stock  [1990], 

While  illustrating  that  spectral  variations  are  theoretically 
possible  at  large  m,  these  arguments  did  not  identify  the 
physical  processes  responsible  for  the  presence  or  absence 
of  an  attenuated  spectrum  in  a  given  region  of  the  atmo¬ 
sphere.  The  first  theory  to  attempt  this  was  developed  by  Kuo 
and  Lue  [1994],  who  proposed  that  the  presence  of  vertical 
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shear  in  the  background  wind  triggered  attenuated  spectra. 
They  argued  that  differential  horizontal  advection  of  vertical 
parcel  oscillations  due  to  the  wave  is  induced  by  the  back¬ 
ground  wind  shear  and  leads  to  reduced  spectral  intensities  in 
the  horizontal-velocity  spectrum.  Zhu  [1994],  on  the  other 
hand,  argued  that  q  <  3  arose  in  the  stratosphere  when 
radiative  damping  exceeded  diffusive  damping  of  the  wave 
spectrum. 

Here  another  theory  is  proposed,  in  which  attenuated 
spectra  arise  due  to  changes  in  background  wind  speeds, 
which  refract  waves  to  smaller  vertical  wavenumbers  and 
produce  accompanying  reductions  in  the  total  shear  variance 
and  vertical  wave-action  flux  of  the  wave  field.  Here  we  shall 
discuss  these  ideas  within  the  context  of  spectral  models  a 
and  b  only,  as  they  are  the  most  obviously  different  from 
each  other  [Hines,  1991a].  In  the  development  it  will  hope¬ 
fully  become  clear  that  this  mechanism  is  sufficiently  general 
and  robust  in  nature  that  it  should  produce  qualitatively 
similar  spectral  responses  within  the  framework  of  model  c 
or  any  other  wave  spectral  theory. 

Indeed,  these  refraction  processes  have  been  encapsu¬ 
lated  mathematically  within  parameterizations  of  gravity- 
wave  processes  for  general  circulation  models  [e.g., 
Lindzen,  1981,  1985;  Fritts  and  Lu,  1993;  Lu  and  Fritts, 
1993].  Thus  for  the  most  part  we  shall  either  apply  or  adapt 
these  formulae  to  illustrate  anticipated  spectral  effects.  The 
major  effort  is  to  present  a  persuasive  case  for  a  connection 
between  theoretical  changes  in  spectral  shape  induced  by 
changes  in  the  background  wind  speed,  and  the  occurrence 
or  otherwise  of  reduced  spectral  densities  at  large  m  in  the 
observational  literature. 

Wind-Shifting  Theory 

Basic  Theory  for  a  Single  Wave 

In  common  with  the  approximations  used  in  developing 
models  a  and  b  [Smith  et  al.,  1987;  Hines,  1991b],  we  shall 
assume  hydrostatic  gravity  waves  in  a  nonrotating  atmo¬ 
sphere.  Such  waves  obey  the  simple  dispersion  relation 
m2  —  k2N2Iw2,  where  k  is  the  horizontal  wavenumber  and 
(o  is  the  intrinsic  frequency.  Given  that  w  =  k(ca  -  U  cos  <p) 
by  definition,  where  ch  is  the  ground-based  horizontal  phase 
speed,  U  is  the  background  horizontal  wind  speed,  and  <p  is 
the  angle  between  »c  and  U,  then 

N 

m  =  ] - fj - r  (D 

\ch  -  U  cos  tp  I 

The  denominator  \ch  —  U  cos  <p\  is  the  intrinsic  horizontal 
phase  speed  of  the  wave,  |c“‘|  =■  |m/k|. 

We  shall  assume  that  horizontal  and  temporal  variations  in 
the  background  atmosphere  are  negligible,  so  that  k  and  ch 
remain  constant  as  a  given  wave  propagates  upward  [e.g., 
Eckermann,  1992].  Since  in  the  absence  of  dissipation  the 
wave  also  conserves  its  Eliassen-Palm  (EP)  flux  p0u'w', 
where  p0  is  the  background  density  and  w'  is  the  vertical 
velocity  perturbation  [Eliassen  and  Palm,  1960],  it  follows 
that  the  horizontal- velocity  variance  u'2  varies  with  height  z 
approximately  as 


where  m0  and  u'0 1  are  the  wavenumber  and  horizontal- 
velocity  variance,  respectively,  at  height  z0,  and  H p  is  the 
density  scale  height  [e.g.,  Lindzen,  1981;  Schoeberl,  1985b]. 
Using  (1)  and  assuming  that  N  is  constant,  (2)  can  be 
reexpressed  as  [Lindzen,  1981] 


where 


11  1  d  _  11 

—  =  —  + - - - U  =  —  + - .  (4) 

Hu  Hp  ch  -  U  cos  <p  dz  Hp  H0 

The  scale  height  for  u'2,  Hu,  quantifies  the  compensating 
interchange  of  variance  between  u'2  and  w' 2  in  order  to 
maintain  both  the  dispersion  relation  (1)  and  a  conserved  EP 
flux.  Similarly,  the  normalized  shear  variance  c r1  =  m2u'2l 
N2,  which  quantifies  the  stability  of  a  wave  (or,  indeed,  of  a 
spectrum  of  waves)  [e.g.,  Desaubies  and  Smith,  1982; 
Hines,  1991a],  varies  with  a  scale  height  Ha  given  by 


11  3  d  _  13 

—  *  —  + - = - U  =  —  + - .  (5) 

H„  Hp  ch  —  U  cos  <p  dz  Hp  H0 

It  is  worth  stressing  at  this  point  that  eyen  though  a 
wind-shear  term  dU/dz  appears  in  (4)  and  (5),  it  plays  no 
direct  physical  role  in  modifying  wave  characteristics 
through  the  momentum  equations,  as  it  can  do  in  other 
situations  [e.g.,  Schoeberl,  1985a;  Hines,  1989].  Rather,  it 
appears  purely  due  to  the  reformulation  of  the  vertical 
variations  of  u'2  as  an  exponential  in  (3).  Note  that  when  (4) 
is  multiplied  by  dz  in  (3),  we  are  left  with  a  term  of  the  form 
A£7  =  £7(z  +  dz)  -  (7(z),  indicating  that  it  is  the  wind 
speed  variations  rather  than  the  shear  itself  that  is  the 
fundamental  process  here.  In  contrast,  the  shear  itself  is 
fundamental  to  the  theory  of  spectral  attenuation  developed 
by  Kuo  and  Lue  [1994]. 

Figure  la  illustrates  how  <p  relates  to  the  phase  velocity  ch 
and  wind  velocity  U  in  the  horizontal  plane.  It  reduces  to  a 
one-dimensional  problem  along  an  axis  QOP  which  is 
aligned  parallel  to  ch,  as  depicted  in  Figure  lb.  Conse¬ 
quently,  we  shall  consider,  without  loss  of  generality,  the 
limiting  case  of  <p  =  0"  or  180°,  so  that 

c'T  =  ch-0,  (6) 

where  C/,  and  U  can  now  be  positive  or  negative  along  the 
QOP  axis  in  Figure  lb.  In  most  cases  that  we  shall  consider, 
U  is  a  zonal  wind,  and  so  QOP  is  aligned  east-west. 

Thus  Hp  can  be  positive  or  negative,  depending  on 
whether  |c  “l|  decreases  or  increases  with  height.  To  monitor 
this,  we  introduce  the  sign  parameter 


with  j 3  defined  such  that  it  matches  the  sign  convention  of  the 
wind  shift  parameter  f3  introduced  later  when  analyzing 
wave  spectra. 

If  p  =  -1,  then  Hp?  <  0.  This  leads,  via  (3)  and  (4), 
to  growth  of  u'2  with  height  which  is  more  gradual  than 
exp  (/  dzIHp),  in  order  to  conserve  EP  flux  as  the  mean 
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wind  varies.  For  p  =  -1,  the  winds  also  dilate  |c™‘|  and  so 
reduce  m  according  to  (1).  This  wavenumber  decrease  is  a 
refraction  effect  induced  by  background  wind  shifts  which 
increase  a>.  Since  we  are  concerned  primarily  with  wave- 
numbers  here,  we  shall  refer  to  j8  =  -1  as  “downshifting.” 

These  effects  are  illustrated  in  Figure  2a  for  the  ch  =  25  m 
s-1  wave  on  the  left,  propagating  in  a  mean  wind  environ¬ 
ment  typical  of  the  winter  lower  stratosphere  at  ~55°S. 

Conversely,  the  other  wave  in  Figure  2a  experiences 
“upshifting”  of  vertical  wavenumbers  m  within  these  back¬ 
ground  winds  (i.e.,  p  =  +1),  since  U  moves  closer  to  this 
wave’s  ch  value  as  it  propagates  upward,  so  that  |c“‘| 
decreases  and  m  increases.  This  upshifting  leads  to  positive 
Hi ;  values  and  thus  variance  growth  with  height  greater  than 

exp  (J  dz/Hp),  as  illustrated.  _ 

Figure  2b  traces  the  evolution  of  u7^  and  m  for  the  ft  =  —  1 
and  ft  =  +1  waves  in  Figure  2a  during  their  propagation 
from  10  to  25  km.  For  the  j3  =  +1  wave,  both  u 77  and  m 
increase,  whereas  both  quantities  decrease  for  the  p  =  —1 
wave.  Effects  such  as  these  were  first  quantified  by  Lindzen 
[1985],  who  investigated  their  influence  on  the  momentum 
budget  of  the  lower  atmosphere. 

Since  background  wind  changes  can  produce  systematic 
changes  in  both  wave  variance  and  wavenumber,  they  can 
influence  the  shape  of  the  vertical  wavenumber  power 
spectrum  of  horizontal  velocities,  Fu(m). 

Application  to  the  Vertical  Wavenumber  Spectrum 
Fritts  and  Lu  [1993]  parameterized  the  effects  of  wind 
shifts  on  a  theoretical  wavenumber  spectrum  based  to  some 
extent  on  a  saturation  model.  However,  the  dynamical 


Figure  2.  (a)  Variation  of  wave  horizontal-velocity  oscilla¬ 
tions  u'  (thin  solid  curve)  with  height  for  two  waves, 
centered  about  their  horizontal  phase  speeds  (shown  with 
dashed  lines)  of  ch  =  25  m  s-1  ^3  =  -1)  and  ch  -  56.9  m 
s-1  (p  =  +1).  The  wind  profile  U  is  the  thick  curve,  and  the 
«'  variations  relative  to  ch  are  shown  on  the  top  x  axis;  exp 
(J  dzHHp)  growth  of  the  initial  velocity  amplitude  is  plotted 
with  dotted  curves,  assuming  an  isothermal  background 
temperature  of  215  K;  (b)  variation  of  u'2  (solid  curve) 
versus  m  for  both  waves  during  propagation  from  10  to  25 
km.  Dotted  curves  show  simple  exp  (J  dz/Hp)  scaling  of  the 
initial  variance  with  height. 


approach  they  used  to  model  the  effects  is  similar  to  the  one 
used  by  Hines  [1991b,  1993b]  to  describe  the  evolution  of  the 
spectrum  under  a  “Doppler-spreading”  model.  The  Dop- 
pler-spreading  effect  is  basically  the  same  as  the  wind- 
induced  refraction  that  we  have  discussed,  except  that  for 
Doppler  spreading.  At/  is  due  to  random  velocity  fluctua¬ 
tions  of  larger-scale  waves  rather  than  a  background  wind. 
This  produces  statistical  spreading  of  wavenumbers  m,  as 
opposed  to  the  systematic  shifts  induced  by  a  discrete 
change  in  the  mean  wind  speed. 

The  analytical  model  of  the  spectrum  assumed  by  Fritts 
and  Lu  [1993]  in  the  absence  of  background  wind  shifts  takes 
the  form 


ch  --■=>  ;  U  cos<p  ••••>■  ;  ch  -U  cos  <p — 

Figure  1.  (a)  Phasor  diagram  showing  the  geometry  of 
ch  —  U  cos  <p  with  varying  azimuth  angle  <p  and  (b^  the 
one-dimensional  simplification  of’tjie  problem  along  QOP. 


a)  Variation  of  Intrinsic  Phase  Speed  withcp 

north 

->  ->  -> 

I  PQ  I  -  I  PO  I  +  I  OQ  I  -  I  ch-  U  cos<P  I 


b)  Intrinsic  Phase  Speed  with  respect  to  QOP 


ch  -U  cosq> 
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Figure  3.  Response  of  F'u(m)  (solid  curve)  to  upshifting 
(dotted  curves,  p_=  +1)  and  downshifting  (dashed  curves, 
p  =  -l),usingAt/  values  of  10and20ms_\lV  =  0.021  rad 
s_1  and  m*  =  2ir/(2.5  km).  The  shaded  region  is  the  excess 
upshifted  spectral  density  which  is  scaled  back  to  the  spec¬ 
tral  saturation  limit  of  aN2m  ~3 . 


m* 


a  N3l4e~zlH>. 


(9) 


This  is  similar  to  the  analytical  Doppler-spread  spectrum 
given  by  Hines  [1991b].  Fritts  and  Lu  [1993]  assumed  that 
the  unshifted  spectrum  F'u(m)  was  invariant  with  altitude  at 
m  »  m*,  so  that  HE  *  2 Hp  =  HJ2  and  F'u(m )  -» 
aN2lm3  at  m  »  m*,  where  a  ~  1/8-1/6  [Fritts  and 
VanZandt,  1993],  in  general  agreement  with  the  unattenu¬ 
ated  observations  cited  earlier.  Under  these  assumptions,  (8) 
is  plotted  as  the  solid  curve  in  Figure  3  using  2trlmif  =2.5 
km  and  N  «  0.021  rad  s_I,  which  are  typical  values  in  the 
lower  stratosphere  [e.g.,  Fritts  and  Chou,  1987;  Tsuda  et  al., 
1991;  Allen  and  Vincent,  1995]. 

To  model  the  effects  of  wind-induced  refraction  on  the 
wave  spectrum,  Fritts  and  Lu  [1993]  considered  the  spectral 
response  to  a  sudden  shift  in  wind  speed  At7,  assuming  no 
change  in  the  altitude  of  observation  and  ascribing  the 
spectral  density  Fi(/n,)  at  some  m,-  to  a  wave  or  waves  of 
wavenumber  m,.  Their  results  for  p  =  +1  and  p  =  -1  are 
shown  with  dotted  and  dashed  curves,  respectively,  in 
Figure  3.  The  form  and  degree  of  spectral  variation  is 
controlled  by  the  wind  shift  parameter 


P  =P 


|Af/|m* 

N 


(10) 


which  is  essentially  the  same  as  the  /3  factor  arising  in  models 
of  Doppler-shifted  frequency  spectra  [e.g.,  Fritts  and  Van¬ 
Zandt,  1987;  Gardner  et  al.,  1993]. 

As  in  Figure  2,  p  =  -1  yields  reduced  wavenumbers 
(downshifting)  and  a  reduction  in  the  area  beneath  the 
spectrum  (i.e.,  reduced  variance).  The  converse  is  true  for 
P  =  + 1,  and  Fritts  and  Lu  [1993]  argued  that  the  increases  in 


total  variance  for  p  =  +1  produced  instabilities  within  the 
wave  field.  For  example,  single-wave  velocity  perturbations 
become  convectively  unstable  when  u'  >  ch  -  U  cos  <p. 
This  occurs  at  a  phase  where  the  —  u'  perturbations  in  Figure 
2a  (which  are  centered  on  ch)  cross  the  0  curve,  as  the  p  = 
+ 1  wave  does  at  z  ~  23.5  km.  Saturation  models  argue  that 
the  amplitudes  of  such  waves  are  attenuated  to  saturated 
levels,  so  that  instabilities  no  longer  occur.  So,  as  waves 
approach  a  critical  level  (|  ch  -  U  cos  <p| — >  0,  m  — *  »),  they 
are  almost  totally  dissipated  («'  — »  0),  although  this  ignores 
various  nonlinear  transmission,  reflection,  and  amplification 
processes  which  can  arise  within  these  regions  [e.g., 
Maslowe,  1986;  Walterscheid  and  Schubert,  1990]. 

On  the  basis  of  these  saturation  arguments,  Fritts  and  Lu 
[1993]  concluded  that  the  increases  in  spectral  density  at 
large  m  for  p  =  +1  (shaded  in  Figure  3)  could  not  be 
sustained  and  so  scaled  their  upshifted  spectra  back  to  the 
aN2/m 3  limit  to  restore  wave  field  stability.  Sato  and 
Yamada  [1994]  applied  saturation  processes  to  a  single 
upshifted  wave  and  derived  model  upshifted  spectra  close  to 
this  limiting  shape. 

Doppler-spreading  models  do  not  attribute  this  limiting 
spectral  shape  to  saturation  but  rather  to  the  characteristic 
statistical  form  that  Doppler-spreading  effects  take  in  the 
Fourier  domain;  but  for  a  spectrum  of  p  =  +1  waves,  wind 
shifts  lead  to  increases  in  both  u7^  and  cr2  over  the  whole 
spectrum,  according  to  (4)  and  (5).  Both  of  these  increases 
trigger  enhanced  Doppler  spreading  of  waves  to  large-m 
values  where  they  are  removed  by  dissipation  [Hines,  1991a, 
b,  1993b].  Once  enough  waves  have  been  removed  so  that  cr2 
is  reduced  to  a  stable  value,  Doppler  spreading  should  return 
the  characteristic  aN2lm  3  spectrum  at  large  m.  Thus  as  long 
as  these  removal  processes  are  sufficiently  rapid,  then  an 
attenuation  of  the  large-m  spectrum  for  p  =  +1  waves, 
which  is  qualitatively  similar  to  that  invoked  by  Fritts  and 
Lu  [1993],  should  also  arise  within  a  Doppler-spreading 
framework. 

Thus  to  summarize,  simple  models  of  refracted  wave 
spectra  predict  reduced  spectral  densities  (a  <  1/10-1/6)  and 
steeper  spectral  indices  (q  >  3)  when  mean-wind  changes 
yield  downshifting  (/?  =  -1),  whereas  conventional 
aA(N)m~3  spectra  are  anticipated  when  mean-wind 
changes  are  small  or  yield  upshifting  0  =  0  or  +1). 

Spectral  Observations  and  the  Background 
Wind  Profile 

We  now  review  various  spectral  measurements  presented 
in  the  literature  and  reanalyze  them  with  respect  to  the 
background  wind  profile  to  test  whether  application  of  these 
ideas  can  account  for  attenuated  and  unattenuated  spectral 
data  at  large  m.  The  survey  of  the  observational  literature 
will  be  selective  rather  than  exhaustive.  The  major  selection 
criteria  for  reanalyzing  a  given  data  set  are  acceptable 
accuracy  (both  instrumental  and  statistical)  and  knowledge 
of  the  background  wind  profiles  accompanying  the  measure¬ 
ments,  to  facilitate  predictions  using  the  wind-shifting 
model.  Such  data  will  range  in  character  from  isolated 
campaign  observations  to  detailed  climatological  studies. 

Lower  Stratosphere 

Radiosonde  and  radar  measurements  at  Kyoto,  Japan.  A 

particularly  thorough  observational  study  of  vertical  wave- 
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number  spectra  in  the  lower  stratosphere  was  made  possible 
by  high-resolution  radar  and  radiosonde  sounding  programs 
conducted  at  a  site  near  Kyoto  (35°N,  136°E)  by  scientists  of 
the  Radio  Atmospheric  Science  Center  of  Kyoto  University. 
The  radar  provided  velocities,  while  the  radiosondes  pro¬ 
vided  temperature  data,  and  initial  intercomparison  of  the 
radar  measurements  of  horizontal-velocity  spectra  and  the 
radiosonde  measurements  of  the  spectrum  of  relative- 
temperature  perturbations  T'/T0  ( T '  =  temperature  pertur¬ 
bation,  T0  —  background  temperature)  revealed  excellent 
agreement  with  theory  [ Fritts  et  al.,  1988].  Similar  analyses 
were  performed  on  more  extensive  data  sets  provided  sub¬ 
sequently  by  these  instruments  at  different  times  of  the  year 
[Tsuda  et  al.,  1989,  1991;  Kuo  and  Lue,  1994]. 

Tsuda  et  al.  [1991]  spectrally  analyzed  over  2  years  of  the 
radiosonde  data.  They  found  that  at  large  m ,  the  spectra  of 
relative-temperature  perturbations  and  wave-perturbed 
Brunt-Vaisala  frequencies  in  the  stratosphere  exhibited  sim¬ 
ilar  seasonal  variations.  Their  model  for  Fj'/To(m)  assumed 
q  —  3  at  large  m,  with  a  polarization  factor  A{N)  = 
0.6N4/g2  and  a  =  1/6,  which  they  compared  with  their 
measured  T'IT0  spectra.  The  results  are  shown  in  Figure  4, 
as  transcribed  from  Figure  12  and  Table  1  of  their  paper. 
While  in  winter  a  values  were  close  to  (although  systemati¬ 
cally  larger  than)  model  predictions  and  q  was  very  close  to 
3,  during  summer  the  spectral  density  at  large  m  was 
consistently  <50%  of  model  predictions,  and  q  =  2. 2-2. 4. 
Kuo  and  Lue  [1994]  analyzed  horizontal-velocity  spectra 
using  a  radar  at  the  same  site  (the  MU  radar)  and  found 
similar  seasonal  differences  in  a.  However,  while  Kuo  and 
Lue  [1994]  found  q  ~  3  in  winter,  their  data  revealed  slightly 
greater  spectral  indices  {q  a  3)  in  summer. 

Figures  5a  and  5b  show  estimates  of  the  mean  zonal  winds 
over  this  site  during  winter  and  summer,  respectively.  The 
stratospheric  values  are  balanced  wind  calculations  using 
climatological  temperature  and  heating  rate  data  from  the 
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Figure  4.  Observed  spectral  indices  q  versus  observed 
power  spectral  densities  at  large  m  (normalized  by  the 
large-m  model  spectral  density  F’TiT<){m))  from  the  radio¬ 
sonde  study  of  Tsuda  et  al.  [1991].  Each  point  is  a  seasonal 
average  in  a  given  year,  and  q  was  obtained  by  fitting  spectra 
over  the  2W(0.5  km)-2W(1.4  km)  wavenumber  interval.  The 
model  prediction  at  large  m  is  asterisked. 
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Figure  5.  Plots  of  zonal  wind  U  versus  height  at  Kyoto 
from  COSPAR  International  Reference  Atmosphere  (CIRA) 
data  (solid  curve)  and  a  schematic  fit  based  on  climatology 
(dotted  line)  (a)  during  winter  and  (b)  summer.  Upward 
pointing  arrows  represent  upward  propagating  waves,  and 
their  phase  speed  ch  is  given  by  the  x  axis  velocity  to  which 
they  point  [e.g.,  Lindzen,  1981].  Waves  for  which  ch  <  U 
have  dark  arrows,  whereas  C/,  >  U  waves  have  lightly 
shaded  arrows.  The  shaded  strip  marks  the  height  interval 
spectrally  analyzed  by  Tsuda  et  al.  [1991]. 


COSPAR  International  Reference  Atmosphere  (CIRA) 
[Marks,  1989],  while  the  tropospheric  profile  is  a  schematic 
continuation  of  these  data  below  the  10-km  CIRA  lower 
boundary  and  is  based  crudely  on  climatology  [Holton,  1992, 
Figure  6.1].  Temperature  profiles  were  constructed  in  the 
same  way  and  were  used  to  convert  the  wind  profiles  from 
log  pressure  to  geometric  heights.  The  height  interval  over 
which  Tsuda  et  al.  [1991]  spectrally  analyzed  their  data  is 
shaded. 

Waves  are  depicted  on  these  diagrams  with  upward  point¬ 
ing  vectors  at  their  ground-based  horizontal  phase  speed  ch, 
as  in  Figure  2a.  As  stated  earlier,  we  assume  that  ch  remains 
constant  with  height,  and  we  choose  its  value  on  the  premise 
that  waves  are  produced  in  the  troposphere,  so  that  ch 
should  be  of  the  order  of  U  values  in  tropospheric  source 
regions.  A  selection  of  waves  for  which  ch  <  U  (e.g., 
mountain  waves,  ch  ~  0)  are  depicted  in  Figure  5a  with 
darkly  arrowed  vectors.  The  wintertime  wind  structure 
allows  most  of  these  waves  to  propagate  through  the  obser¬ 
vational  height  range,  as  is  well  known  [e.g.,  Lindzen,  1985; 
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Schoeberl,  1985a],  The  background  wind  structure  produces 
downshifting  0  =  -1)  of  these  waves  in  the  troposphere, 
since  dU/dz  >  0  and  so  U  moves  away  from  ch  with 
increasing  altitude,  whereupon  e“‘  increases  and  m  de¬ 
creases.  Conversely,  in  the  stratosphere,  dU/dz  <  0,  lead¬ 
ing  to  upshifting  0  =  +1)  since  U  moves  nearer  to  ch  as  the 
wave  propagates  upward  and  m  increases,  as  in  Figure  2a. 
Indeed,  Sato  and  Yamada  [1994]  presented  direct  evidence 
of  upshifted  vertical  wavenumbers  in  radiosonde  and  MU 
radar  data  from  these  heights  during  October,  when  the 
background-wind  profile  is  similar  to  Figure  5a. 

Waves  for  which  ch  >  U  are  shown  with  lightly  shaded 
arrows  in  Figure  5a.  Lower-tropospheric  waves  for  which  10 
m  s~‘  ^  ch  s50ms"’  cannot  enter  the  stratosphere  since 
they  encounter  critical  levels.  However,  observations  sug¬ 
gest  that  gravity  waves  are  also  produced  near  the  jet  stream 
peak  at  ~10  km  [e.g.,  Hooke  and  Hardy,  1975;  Uccellini  and 
Koch,  1987;  Yamanaka  et  al.,  1989;  Fritts  and  Nastrom, 
1992;  Sato,  1994],  and  these  waves  can  enter  the  strato¬ 
sphere  where  they  experience  downshifting  0  =  - 1). 

Observations  and  modeling  indicate  that  over  undular  land 
masses,  topographic  forcing  of  gravity  waves  is  more  intense 
than  production  by  other  sources  [e.g.,  Rind  et  al.,  1988; 
Nastrom  and  Fritts,  1992;  Fritts  and  Nastrom,  1992;  Bac- 
meister,  1993].  On  analyzing  4  years  of  MU  radar  measure¬ 
ments  at  heights  —  16  km,  Murayama  et  al.  [1994]  found  u'w' 
<  0  in  winter,  consistent  with  ch  <  U  waves.  Similar 
wintertime  directions  were  inferred  by  Kitamura  and  Hirota 
[1989]  from  stratospheric  radiosonde  profiles  over  Japan. 
Sato  [1994]  presented  strong  evidence  that  the  bulk  of  the 
velocity  variance  observed  by  the  MU  radar  at  —18-22  km 
during  winter  was  due  to  mountain  waves  with  ch  values 
near  zero.  This  indicates  that  mountain  waves  (dark  arrows) 
are  more  energetic  than  waves  emanating  from  the  jet  stream 
(lightly  shaded  arrows)  during  winter  at  this  site.  However, 
even  if  they  had  similar  EP  fluxes,  Figure  5a  reveals  that  the 
mountain  waves  undergo  upshifting  0  =  +1)  in  the  strato¬ 
sphere  while  waves  from  near  the  jet  stream  are  downshifted 
0  =  -1).  This  differential  refraction  yields  amplification  of 
the  mountain-wave  velocity  variance  relative  to  the  variance 
of  jet-stream-associated  waves,  as  depicted  in  Figure  2a. 

Thus  the  stratospheric  wave  variance  during  winter  should 
be  dominated  by  waves  for  which  ch  —  0.  Since  the  decrease 
of  the  zonal  wind  with  height  produces  upshifting  effects 
0  =  +1)  on  such  a  wave  field,  this  will  give  rise  to  wave 
dissipation  and  a  forcing  of  the  spectrum  to  the  aA(N)m  -3 
limit  at  large  m,  according  to  the  wind-shifting  model.  The 
data  of  Tsuda  et  al.  [1991]  support  these  predictions  (see 
Figure  4). 

Figure  5b  repeats  this  analysis  for  the  summertime  wind 
profile.  In  this  case,  the  ch  <  U  waves,  such  as  mountain 
waves,  cannot  propagate  into  the  height  range  analyzed  by 
Tsuda  et  al.  [1991],  due  to  critical-level  encounters  lower 
down.  Murayama  et  al.  [1994]  found  u'w'  ~  0  in  summer¬ 
time  MU  radar  observations  at  —16  km,  a  height  at  which 
mountain  waves  can  still  exist  in  the  data  (see  Figure  5b). 
Since  u'w'  <  0  for  mountain  waves,  the  near-zero  u'w' 
values  measured  by  Murayama  et  al.  [1994]  suggest  that 
ch>  U  waves  also  exist  at  —16  km  with  approximately 
equal  and  opposite  EP  fluxes.  On  propagating  beyond  z  =  16 
km,  Figure  5b  shows  that  many  waves  for  which  u'w'  <  0 
reach  critical  levels,  leaving  mainly  ch  >  U  (u’w'  >  0) 
waves  at  the  heights  analyzed  by  Tsuda  et  al.  [1991]. 


[Tsuda  etal.,  1994a,  1994 b\ 
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Figure  6.  (a)  Mean  zonal  and  meridional  velocity  spectra 
from  23  to  36  km  over  Watukosek,  Indonesia  (7.6°S, 
112.7°E),  from  February  27  to  March  22,  1990.  The  model 
spectrum  of  Smith  et  al.  [1987]  at  large  m  is  also  plotted,  (b) 
Mean  zonal  wind  as  a  function  of  height  from  the  radiosonde 
data  during  this  period.  As  in  Figure  5,  the  arrows  represent 
waves,  while  the  shaded  region  marks  the  23-  to  36-km  band 
spanned  by  the  data  in  Figure  6a.  The  data  were  transcribed 
from  Tsuda  et  al.  [1994a,  b]. 


Thus  in  summer  the  stratospheric  wave  field  should  be 
dominated  by  the  eastward  phase-speed  waves  from  the  jet 
stream,  for  which  j3  =  - 1  since  U  moves  away  from  these  ch 
values  as  they  propagate  through  the  shaded  height  range.  In 
these  circumstances  the  wind-shifting  model  predicts  re¬ 
duced  spectral  densities  at  large  m,  as  observed  by  Tsuda  et 
al.  [1991]  and  Kuo  and  Lue  [1994].  However,  the  Fritts  and 
Lu  [1993]  model  predicts  spectral  indices  q  >  3,  which  agree 
with  the  radar  results  of  Kuo  and  Lue  [1994]  but  not  the 
radiosonde  findings  of  Tsuda  et  al.  [1991], 

Equatorial  measurements  over  Indonesia.  Tsuda  et  al. 
[1994a,  b]  analyzed  radiosonde  measurements  of  horizontal 
velocities  and  temperatures  from  a  month-long  observational 
campaign  at  Watukosek,  Indonesia  (7.6°S,  112.7°E).  Tsuda 
et  al.  [1994b]  noted  that  stratospheric  horizontal-velocity 
spectra  at  large  m  were  only  10-25%  of  model  predictions 
and  that  q  >  3.  These  spectra  are  reproduced  in  Figure  6a. 

Figure  6b  plots  the  mean  background  zonal  winds,  derived 
from  all  the  measurements  during  the  campaign  [Tsuda  et 
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(a)  Tsuda  et  al.  [1991]:  WINTER 


(b)  Tsuda  etal.  [1991]:  SUMMER 
401 . . i . . \ 


I 


_  1  10 
7*J/(model) 


(c)  Tsuda  etal.  [1994a,  1994b] 


Figure  7.  Plots  of  T 2  in  the  150-  to  900-m  range  of  vertical  scales,  normalized  by  model-predicted 
variances,  (a)  and  (b)  Seasonal-mean  profiles  during  winter  and  summer,  respectively,  from  the 
observational  study  of  Tsuda  et  al.  [1991].  (c)  Mean  results  from  the  campaign  study  of  Tsuda  et  al. 
[1994b],  The  shaded  strips  show  the  height  intervals  over  which  data  were  spectrally  analyzed,  yielding  the 
results  in  Figures  4  and  6a.  The  dashed  curve  shows  model  predictions. 


al.,  1994a,  Figure  2].  Following  the  presentation  in  Figure  5, 
the  height  interval  spectrally  analyzed  by  Tsuda  et  al. 
[1994b]  is  shaded,  and  waves  for  which  ch  exceeds  (is  less 
than)  €  are  plotted  with  dark  (light  shaded)  vector  arrows. 
The  wind-shifting  environment  is  qualitatively  similar  to  that 
in  Figure  5b,  where  stationary  and  westward  ch  waves  reach 
critical  levels,  so  that  mainly  eastward  phase  speed  waves 
reach  the  stratosphere. 

In  the  23-  to  36-km  height  range,  U  increases  from  ~23  to 
27  km,  then  decreases  from  ~27  to  36  km.  For  eastward  ck 
waves  this  yields  /3  =  + 1  at  lower  heights  and  (3  =  - 1  at  the 
upper  heights,  and  so  /3  =  - 1  dominates  over  most  of  the  23- 
to  36-km  interval.  Therefore  we  anticipate  a  net  downshift¬ 
ing  of  the  spectrum  (i.e.,  (3  =  —1),  for  which  the  wind- 
shifting  model  predicts  reduced  spectral  densities  and  q  >  3 
at  large  m,  as  observed  in  Figure  6a.  However,  strong 
planetary  wave  activity  in  these  data  [Tsuda  et  al.,  1994a] 
may  also  produce  time-varying  modulation  of  the  wave 
transmission  and  wind  shifting  and  hence  time-varying  mod¬ 
ifications  of  the  spectral  shape  [e.g.,  Lu  and  Fritts,  1993]. 

Allen  and  Vincent  [1995]  argued  that  temperature  spectra 
measured  by  radiosondes  in  the  stratosphere  can  be  dis¬ 
torted  by  slow  sensor  response  times,  whereas  variances 
tended  to  be  less  affected.  So,  as  a  further  check  of  these 
radiosonde  data,  Figure  7  shows  height  profiles  of  the 
temperature  variances  at  vertical  scales  between  150  and  900 
m.  These  variances  were  normalized  by  the  model  variance, 
obtained  by  integrating  FlT(m)  over  the  same  wavenumber 
interval  [Tsuda  et  al.,  1991,  1994b].  The  dashed  line  shows 
model  predictions.  We  note  that  the  height  variations  in  the 
shaded  regions  of  Figures  7b  and  7c  compare  very  well  with 
the  refractive  attenuation  of  variance  predicted  using  the 
wind  profiles  in  Figures  5b  and  6,  respectively.  Variances  in 
Figure  7b  attenuate  above  —18  km,  a  level  at  which  moun¬ 
tain  waves  are  removed  by  critical  levels,  leaving  mainly 
downshifted  jet  stream  waves  above  this  level  (see  Figure 
5b).  Variances  in  Figure  7c  only  fall  below  model  predictions 
above  —28  km,  the  same  height  above  which  downshifting 
was  predicted  in  Figure  6.  The  winter  profiles  in  Figure  7a 
are  near  model  values  throughout  most  of  the  shaded  region. 


again  consistent  with  earlier  predictions  (see  Figure  5a). 
However,  values  attenuate  somewhat  at  the  uppermost 
heights.  This  could  indicate  removal  of  ck  <  U  waves  at 
these  heights,  yielding  downshifting  of  ck  >  U  waves,  since 
zonal  wind  profiles  measured  by  Kitamura  and  Hirota  [1989] 
and  Murayama  et  al.  [1994]  over  Japan  often  reveal  U  near 
zero  at  heights  —23-25  km. 

Radiosonde  measurements  in  the  southern  hemisphere. 
Allen  and  Vincent  [1995]  undertook  a  detailed  statistical 
study  of  lower-stratospheric  spectra  using  high-resolution 
radiosonde  temperature  measurements  from  around  Austra¬ 
lia  and  Antarctica.  Unlike  the  results  of  Tsuda  et  al.  [1991] 
from  a  near-conjugate  location  in  the  northern  hemisphere, 
Allen  and  Vincent  [1995]  noted  little  seasonal  variability  in 
their  spectra.  The  only  location  where  spectral  densities  at 
large  m  departed  significantly  from  the  aA{N)m  ~3  limit  was 
at  Davis,  Antarctica  (69°S,  78°E),  where  reduced  spectral 
densities  and  q  <  3  were  observed.  All  of  their  spectra  were 
corrected  for  sensor-response  distortion. 

Figure  8  shows  CIRA  zonal  wind  profiles  during  summer 
and  winter  over  two  of  their  radiosonde  sites:  Davis,  Ant¬ 
arctica  (69°S,  78°E)  and  Woomera  (31°S,  136°E),  a  station 
which  lies  roughly  in  the  middle  of  the  various  Australian 
radiosonde-launching  sites.  The  height  interval  spectrally 
analyzed  by  Allen  and  Vincent  [1995]  is  shaded.  Figures  8b 
and  8d  reveal  that  in  the  winter  stratosphere,  dU/dz  differs 
in  sign  between  Davis  and  Woomera.  Thus  during  winter, 
downshifting  (j8  =  —1)  dominates  at  Davis,  while  upshifting 
dominates  at  Woomera  if  ch  =  0  waves  exist  (due  to  the 
amplification  of  upshifted  variance  relative  to  the  down¬ 
shifted  variance).  This  is  consistent  with  the  attenuated 
spectra  observed  at  Davis  and  unattenuated  spectra  ob¬ 
served  at  the  Australian  stations  by  Allen  and  Vincent 
[1995].  Spectral  indices  q  were  less  than  3  at  Davis,  whereas 
the  wind-shifting  model  predicts  q  >  3. 

Except  for  Davis,  all  data  analyzed  by  Allen  and  Vincent 
[1995]  came  from  sites  equatorward  of  45°S,  and  all  the 
spectra  from  these  sites  tended  to  aA(N)m  ”3  at  large  m,  in 
contrast  to  the  attenuated  Davis  spectra.  Figure  9  shows  a 
latitude-height  cross  section  of  U  at  135°E  using  CIRA  July 
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January  (Summer) 


a)  Davis 


July  (Winter) 
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Figure  8.  Climatological  CIRA  zonal  winds  during  (a)  and 
(c)  summer  and  (b)  and  (d)  winter  over  Davis,  Antarctica, 
and  Woomera,  Australia.  As  before,  arrows  represent 
waves,  while  shading  marks  the  regions  spectrally  analyzed 
by  Allen  and  Vincent  [1995], 


(winter)  data.  The  diagram  shows  that  U  a  0  throughout 
most  of  the  lower  stratosphere  and  that  dU/dz  <  0  equator- 
ward  of  45°S,  whereas  dU/dz  >  0  poleward  of  45°S. 
According  to  the  wind-shifting  model,  this  will  yield  down¬ 
shifted  (attenuated)  spectra  poleward  of  45°S  (Davis)  and 
upshifted  (unattenuated)  spectra  equatorward  of  45°S  during 
winter,  as  observed  by  Allen  and  Vincent  [1995].  This 
assumes  that  the  winter  ch  distributions  do  not  vary  too 
much  with  latitude. 

Figures  8a  and  8c  show  the  wind  profiles  during  summer, 
and  in  these  cases  wind-shifting  predictions  are  less  straight¬ 
forward.  At  Woomera,  ch  =  0  waves  enter  the  shaded 
height  interval  but  probably  cannot  propagate  much  beyond 
20  km.  Hence  a  downshifted  spectrum  dominated  by  the 
lightly  shaded  arrows  might  be  expected,  at  least  above  —20 
km,  whereas  Allen  and  Vincent  [1995]  found  a  seasonally 
invariant  shape  near  cxA(N)m~ 3  at  Woomera.  The  summer¬ 
time  winds  at  Davis  are  quite  weak,  and  so  one  might 
anticipate  a  spectrum  close  to  F‘u(m),  yet  Allen  and  Vincent 
found  attenuated  spectra  which  were  similar  in  shape  to  the 
winter  spectra. 

Analysis  of  these  summertime  spectra  using  the  wind- 
shifting  model  is  made  difficult  by  a  lack  of  knowledge  of 
source  climatologies  (and  hence  ch  values)  in  the  southern 


hemisphere.  Topographic  forcing  is  thought  to  be  less  in¬ 
tense  in  the  southern  hemisphere,  although  mountainous 
regions  in  Australia  and  Antarctica  do  produce  some  strato¬ 
spheric  wave  activity  [e.g.,  Bacmeister  et  al.,  1990;  Bac- 
meister,  1993].  Thus  wave  production  by  meteorological 
phenomena  in  the  southern  hemisphere  [e.g.,  King  et  al., 
1987;  Eckermann  and  Vincent,  1993;  Pfister  et  al.,  1993; 
Marks  and  Revell,  1994;  Smith  et  al.,  1995]  may  be  more 
important.  Until  more  is  known,  our  analysis  of  these 
summertime  spectral  findings  remains  uncertain. 

Upper  Stratosphere  and  Lower  Mesosphere 

Lidar  measurements  in  Wales  and  France.  Climatological 
data  on  wave  spectra  at  heights  —30-60  km  have  been 
derived  from  Rayleigh  lidar  observations  at  Aberystwyth, 
Wales  (52°N,  4°W),  and  at  two  French  sites,  Biscarosse 
(44°N,  1°W)  and  Haute  Provence  (44°N,  6°E).  Case  studies 
of  typical  spectral  shapes  were  presented  by  Mitchell  et  al. 
[1994]  and  Wilson  et  al.  [1990,  1991a],  and  a  climatological 
picture  of  the  space-time  variability  of  spectra  from  several 
years  of  observations  with  the  Welsh  and  French  lidars  was 
given  by  Marsh  et  al.  [1991]  and  Wilson  et  al.  [1991b], 
respectively. 

Figure  10  plots  CIRA  zonal  winds  at  42.5°N,  0°E  using 
CIRA-balanced  winds.  The  three  height  intervals  spectrally 
analyzed  by  Wilson  et  al.  [1990,  1991b]  are  shaded  sepa¬ 
rately.  The  diagram  shows  that  intensifying  stratospheric 
winds  produce  p  =  —  1  waves  in  the  upper  stratosphere  due 
to  critical-level  removal  of  j8  =  +1  waves,  as  has  been 
verified  observationally  [e.g.,  Hamilton,  1991 ;  Eckermann  et 
al.,  1995]. 

According  to  the  wind-shifting  model,  this  should  yield 
strongly  downshifted  spectra  and  attenuated  spectral  densi¬ 
ties  at  large  m  in  both  summer  and  winter.  The  data  in  the 
height  range  30-45  km  revealed  spectral  densities  at  large  m 
an  order  of  magnitude  lower  than  F'(m)  in  summer  and 
winter  [Wilson  et  al.,  1991b],  in  line  with  these  predictions. 
At  large  m,  spectral  indices  q  were  —3  in  both  seasons, 
although  at  smaller  m,  the  q  values  were  smaller  and 
exhibited  seasonal  variations  [Marsh  et  al.,  1991;  Wilson  et 
al.,  1991b], 

The  reductions  in  spectral  density  due  to  wind  shifting 
must  compete  against  exp  (f  dzlHp)  growth  of  u'w'  for  all 
the  waves  in  the  spectrum.  While  observed  spectral  densities 
at  large  m  were  larger  at  45-60  km  than  at  30-45  km  [Wilson 
et  al.,  1991b],  they  were  still  below  the  aA{N)m ~3  limit. 
This  can  occur  because  downshifting  is  still  significant  (i.e., 
P  =  - 1  and  |/?|  »  0)  in  this  height  interval,  as  evident  from 
Figure  10.  Eventually,  exp  (J  dzlHp)  growth  will  produce 
gravity  wave  instability  and  dissipation,  and  the  observed 
reversal  of  dU/dz  in  the  upper  height  range  60-75  km  in 
Figure  10  is  the  well-known  sign  that  gravity  wave  dissipa¬ 
tion  is  occurring  at  these  levels  and  above  [e.g.,  Lindzen, 
1981;  Holton,  1983].  Dissipation  implies  a  spectrum 
aA(N)m -3  at  large  m  according  to  the  saturation  model 
[e.g.,  Smith  et  al.,  1987;  Fritts  and  Lu,  1993],  while  for  the 
Doppler-spreading  model,  it  implies  that  wavenumbers  are 
spread  to  dissipative  scales  at  large  m  and  that  the  spectrum 
relaxes  to  the  aA(N)m~ 3  limit  [Hines,  1991b]. 

The  reversal  in  dU/dz  due  to  dissipation  of  certain  waves 
forces  upshifting  (j 3  =  +1)  upon  the  full  spectrum  of  waves 
and  thus  according  to  the  wind-shifting  model  will  continu¬ 
ally  force  the  spectrum  to  the  aA{N)m  ~3  limit  at  lower 
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Figure  9.  Latitude-height  cross  section  of  the  CIRA  climatological  zonal  winds  in  the  lower  stratosphere 
at  135°E  ±  12.25°  in  the  southern  hemisphere  during  July  (winter).  Darkest  shading  indicates  strong 
eastward  winds  and  the  lightest  shading  indicates  winds  near  zero  (see  contour  key). 


mesospheric  heights.  Wilson  et  al.  [1991b]  found  that  their 
spectra  at  60-75  km  attained  this  limiting  form  in  both  winter 
and  summer,  in  agreement  with  these  wind-shifting  predic¬ 
tions. 

Case  studies  from  other  sites.  Shibata  et  al.  [1988]  ana¬ 
lyzed  spectra  on  12  nights  throughout  a  year  at  Fukuoka, 
Japan  (34°N,  130°E).  Nine  of  the  nights  analyzed  were  from 
September  to  December,  and  they  presented  mean  T'IT0 
spectra  within  the  upper  stratosphere  (30-49.2  km)  and 
lower  mesosphere  (45-64.2  km).  Tsuda  et  al.  [1992]  studied 
horizontal- velocity  spectra  in  the  height  intervals  22-32  km, 
33-43  km,  and  47-56  km  from  rocketsonde  experiments  over 
Japan  during  January-February  1990. 

The  winds  over  this  site  during  winter  are  similar  to  Figure 
10a,  except  that  the  tropospheric  jet  is  stronger,  winds 
decrease  to  near  zero  at  25-30  km  (see  Figure  5a),  then 
increase  to  a  peak  —80  ms-1  at  65  km  [Marks,  1989]  before 
decreasing  with  height  thereafter.  However,  one  should  also 
bear  in  mind  that  stratospheric  warmings  in  the  winter 
northern  hemisphere  can  drastically  alter  this  climatological 
wind  structure  [e.g.,  Labitzke,  1981],  So,  using  the  winds  in 
Figure  10a  as  a  guide,  we  anticipate  attenuated  spectra 
according  to  the  wind-shifting  model,  particularly  in  the  30- 
to  50-km  range.  Both  Shibata  et  al.  [1988]  and  Tsuda  et  al. 
[1992]  observed  attenuated  spectral  densities  in  this  height 
range.  At  large  m,  Shibata  et  al.  [1988]  found  q  —  3, 
whereas  the  data  of  Tsuda  et  al.  [1992]  indicated  q  >  3.  At 
the  upper  height  interval,  spectra  from  the  lidar  (45-64.2  km) 
and  rocketsonde  (47-56  km)  were  closer  to  the  limiting  form 
otA(N)m  _3 ,  consistent  with  the  reversal  of  the  wind  shear  at 
—65  km,  and  the  ensuing  onset  of  upshifting  (/3  =  +1),  wave 
dissipation,  and  limiting  spectral  forms  at  these  heights 
according  to  the  wind-shifting  model. 

We  finish  here  with  some  spectral  data  from  Arecibo 
(18°N,  67°W)  that  are  less  easy  to  interpret.  Spectra  of 


relative-density  perturbations  p'/p0  at  20-45  km,  measured 
by  a  Rayleigh  lidar  at  Arecibo  during  January,  March,  and 
April  1989,  were  reported  by  Beatty  et  al.  [1992]  and  Senft  et 
al.  [1993].  They  found  attenuated  spectra  at  large  m  and 
spectral  indices  q  —  3  or  slightly  less.  Tsuda  et  al.  [1989] 
computed  horizontal-velocity  spectra  from  Arecibo  UHF 
radar  dati  at  —13-22  km  during  June  1983  and  found  spectra 
close  to  the  aA(N)m  ~3  limit. 

The  climatological  wind  profiles  for  these  months  are 
shown  in  Figure  11.  In  both  cases,  waves  for  which  ch  <  U 
and  ch  >  U  enter  the  spectral  height  interval.  There  is 
evidence  that  mountain  waves  generated  around  Arecibo  are 
observed  in  the  radar  data  [Hines,  1989;  Cornish  and 
Larsen,  1989].  Assuming  that  they  are  more  energetic  than 
jet-stream-associated  waves  [e.g.,  Fritts  and  Nastrom, 
1992],  then  downshifting  dominates  in  Figure  11a  and  up- 
shifting  dominates  in  Figure  lib,  yielding  spectral  predic¬ 
tions  consistent  with  the  observations.  However,  this  pre¬ 
sumes  that  the  upshifting  of  any  ch  >  U  waves  at  upper 
heights  in  Figure  11a  is  not  energetically  significant,  other¬ 
wise  unattenuated  spectra  might  be  expected  in  both  cases. 

Mesosphere 

In  the  mesosphere,  where  the  wind  shear  reverses  (see 
Figure  10),  Senft  et  al.  [1993]  found  spectra  nearer  the 
aA(N)m~ 3  limit,  consistent  with  wind-shifting  predictions 
(i.e.,  a  change  from  stratospheric  downshifting  to  mesos¬ 
pheric  upshifting).  Hostetler  and  Gardner  [1994]  also  ob¬ 
served  attenuated  stratospheric  spectra  and  unattenuated 
mesospheric  spectra  from  airborne  lidar  observations  near 
Hawaii.  We  shall  now  analyze  more  extensive  databases  of 
mesospheric  spectra  to  provide  better  assessments  of  their 
relationship  to  wind-shifting  predictions. 

Most  observations  of  mesospheric  spectra  to  date  have 
revealed  shapes  in  agreement  with  the  aA(N)m  _3  limit  at 
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Figure  10.  Climatological  CIRA  zonal  winds  over  southern 
France  (42.5°N  ±  2.5°,  0°E  ±  12.25°),  using  the  same 
phenomenological  representation  as  Figure  5.  The  three 
shaded  strips  correspond  to  the  three  regions  spectrally 
analyzed  by  Wilson  et  al.  [1991b]. 


large  m  [e.g.,  Smith  et  al.,  1987;  Tsuda  et  al.,  1989;  Wu  and 
Widdel,  1989,  1991].  Since  gravity  wave  dissipation  in  the 
mesosphere  forces  background  wind  speeds  to  decrease,  and 
critical-level  filtering  removes  ch  >  U  waves  in  the  strato¬ 
sphere,  then  mesospheric  gravity  waves  always  tend  to  be 
upshifted  (/3  =  +1;  see  Figure  10).  In  such  cases  the 
wind-shifting  model  predicts  the  limiting  form  aA{N)m  ~3  at 
large  m,  in  accord  with  observations. 

Yet  certain  situations  may  give  rise  to  wind  profiles  in 
which  upshifting  does  not  always  occur  in  the  mesosphere 
[e.g.  ,  Lu  and  Fritts,  1993],  According  to  the  wind-shifting 
model,  observations  at  times  of  downshifting  can  produce 
attenuated  spectra,  even  at  these  heights.  However,  because 
of  its  height  dependence  (equation  (9)),  m*  is  smaller  at 
these  heights,  so  that  wind  shifts  A  £7  need  to  be  correspond¬ 
ingly  larger  to  produce  a  given  degree  of  spectral  downshift¬ 
ing  p,  according  to  (10).  Spectral  measurements  at  times  of 
{significantly  negative  p  in  the  mesosphere  would  provide  a 
powerful  test  of  the  wind-shifting  theory. 

A  climatology  of  lidar  spectra  at  large  m  in  the  mesosphere 
over  Urbana  (40°N,  88°W)  was  presented  by  Senft  and 
Gardner  [1991].  They  noted  greater  variability  in  q  and  a 
than  allowed  for  within  saturation  models.  This  variability, 
while  large,  tended  to  be  quasi  random.  Simultaneous  wind 


data  were  not  available  for  direct  wind-shifting  predictions. 
Climatological  wind  patterns  measured  over  Urbana  [ Franke 
and  Thorsen,  1993]  suggest  upshifted  unattenuated  spectra. 

According  to  wind-shifting  theory,  this  transient  variabil¬ 
ity  suggests  that  the  major  wind-shifting  influences  occur 
over  intraseasonal  timescales  (i.e.,  =£1-2  months).  Franke 
and  Thorsen  [1993]  noted  strong  tidal  activity  at  these 
heights,  which  Lu  and  Fritts  [1993]  showed  can  have  strong 
modulating  interactions  with  the  gravity  wave  spectrum. 
Using  m*  =  2ir/(13  km)  [Senft  and  Gardner,  1991],  N  ~ 
0.02  rad  s-1  [Senft  et  al.,  1994],  and  typical  tidal  variations 
about  a  mean  wind  profile,  A  £7,  of  ±20  ms-1  [Franke  and 
Thorsen,  1993],  one  obtains  j8  variations  of  ~±0.5,  which 
may  account  for  some  of  the  quasi-random  spectral  variabil¬ 
ity  observed  by  Senft  and  Gardner  [1991].  Thayaparan  et  al. 
[1995]  presented  direct  evidence  of  tidally  modulated  gravity 
wave  variances  using  mesospheric  wind  data  from  a  nearby 
site. 

Wu  and  Widdel  [1989,  1991]  observed  aA(N)m~3  spectra 
at  large  m  on  spectrally  analyzing  high-resolution  foil-chaff 
data  from  —80-95  km  at  Andpya  (69°N,  16°E)  during  July  15, 
1987  (summer).  Background  winds  measured  during  these 
observations  [Lubken  et  al.,  1990]  are  shown  in  Figure  12b. 
They  reveal  filtering  characteristics  and  mesospheric  wind 
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Figure  11.  Climatological  zonal  winds  over  the  region 
17.5°N  ±  2.5°,  67.5°W  ±  12.25°  (near  Arecibo)  for  (a)  March 
and  (b)  June.  The  height  intervals  spectrally  analyzed  by 
Senft  et  al.  [1993]  and  Tsuda  et  al.  [1989]  are  shaded  in 
Figures  11a  and  lib,  respectively.  Arrows  have  the  same 
meaning  as  in  Figure  5. 
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Mean  Winds  over  Andaya 
(a)  Winter  [Meyer  et  al.,  1987] 


(b)  Summer  [Lubken  etaL,  1990] 


Figure  12.  (a)  Zonal  (solid)  and  meridional  (dotted)  winds 
over  And0ya  (69°N,  16°E)  for  January  1984  [Meyer  et  al., 
1987].  The  height  interval  spectrally  analyzed  by  Wu  and 
Widdel  [1990]  is  shaded,  (b)  Zonal  winds  at  the  same  location 
during  July  1987  [. Lubken  et  al.,  1990].  The  height  interval 
spectrally  analyzed  by  Wu  and  Widdel  [1989,  1991]  is 
shaded.  CIRA  climatological  zonal  winds  are  shown  with 
dashed  curves  in  each  plot.  Arrows  have  the  same  meaning 
as  in  Figure  5. 


patterns  which  should  yield  upshifted  (unattenuated)  spec¬ 
tra,  in  line  with  Wu  and  Widdel’ s  findings. 

However,  their  analysis  of  data  from  similar  experiments 
during  January-February  1984  (winter)  at  heights  —72-88  km 
[Wm  and  Widdel,  1990]  revealed  a  mean  q  —  2.5  (as  fitted 
over  the  vertical-scale  range  0.1-2.1  km),  although  spectral 
shapes  and  intensities  varied  considerably  among  the  mea¬ 
surements.  The  mean  January  background  winds  measured 
during  these  experiments  [Meyer  et  al.,  1987]  are  shown  in 
Figure  12a  along  with  a  climatological  CIRA  profile  [Marks, 
1989].  Zonal  winds  decreased  to  —79  km,  then  increased 
with  height  thereafter,  which  should  produce  downshifting 
above  79  km.  Furthermore,  mean  meridional  winds  in¬ 
creased  from  near  zero  at  70  km  to  —20  ms-1  at  80  km, 
which  should  accentuate  a  downshifting  trend  by  yielding 
increasing  northwestward  winds  acting  on  a  filtered  mesos¬ 
pheric  wave  field  containing  mainly  south-westward  ch 
values. 

The  form  of  the  mean  wind  profiles  [Meyer  et  al.,  1987] 


raises  the  possibility  that  the  observation  of  q  —  2.5  might 
reflect  attenuated  (downshifted)  mesospheric  spectra  at 
these  heights  and  times,  in  the  same  way  that  smaller  q 
values  often  accompany  spectral  attenuation  in  the  strato¬ 
sphere  [e.g.,  Tsuda  et  al.,  1991].  Furthermore,  wind  speed 
increases  above  —79  km  themselves  may  indicate  reduced 
gravity  wave  dissipation  (and  hence  unsaturated  down¬ 
shifted  spectra),  since  wind  speed  decreases  with  height  are 
driven  by  gravity  wave  dissipation  [e.g.,  Holton,  1983]. 
Cause  and  effect  become  somewhat  blurred  in  such  a  pic¬ 
ture,  but  it  suggests  an  intimate  association  between  upshift- 
ing  (downshifting)  and  wave  dissipation  (stability)  in  the 
mesosphere. 

However,  while  Wu  and  Widdel  [1990]  did  not  compare 
their  spectra  with  the  limiting  a A{N)m  -3  shape,  inspection 
suggests  that  their  mean  spectra  do  not  lie  below  this  model 
curve.  Since  there  is  greater  wave  activity  in  these  winter 
data  than  during  summer  [e.g.,  von  Zahn  and  Meyer,  1989], 
it  seems  likely  that  q  <  3  may  often  arise  here  due  to  the 
existence  of  inertial-range  turbulence  ( q  —  5/3)  at  scales 
s 300-600  m  [e.g.,  Hillert  et  al.,  1994].  Nevertheless,  the 
spectral  shapes  and  intensities  measured  by  Wu  and  Widdel 
[1990]  varied  appreciably  from  measurement  to  measure¬ 
ment.  Background  winds  also  varied  considerably  from 
profile  to  profile  due  to  stratospheric  warmings  and  tides 
[Meyer  et  al.,  1987;  Rottger  and  Meyer,  1987],  thus  yielding 
time-varying  wind  shifts  which  may  be  responsible  for  some  of 
this  transient  spectral  variability  [e.g.,  Lu  and  Fritts,  1993]. 

The  interpretation  of  these  mesospheric  data,  while  admit¬ 
tedly  speculative  and  inconclusive,  at  least  serves  to  high¬ 
light  the  potential  for  attenuated  mesospheric  spectra  ac¬ 
cording  to  the  wind-shifting  model. 

Troposphere 

As  in  the  mesosphere,  spectral  analyses  of  tropospheric 
data  generally  reveal  an  unattenuated  aA(N)m  ~3  shape  at 
large  m  [e.g.,  Tsuda  et  al.,  1991],  Comparison  of  these 
observations  with  wind-shifting  predictions  is  difficult  be¬ 
cause  this  is  the  region  where  most  waves  are  produced,  and 
so  little  critical-level  filtering  of  ch  values  has  yet  taken  place 
to  enable  us  to  identify  whether  upshifting  or  downshifting  is 
preferred.  In  many  situations,  both  processes  probably 
occur,  in  which  case  the  upshifted  waves  will  tend  to  be 
enhanced  in  variance  relative  to  the  downshifted  waves  (see 
Figure  2).  A  spectrum  whose  variance  is  dominated  by 
upshifted  waves  should  produce  an  aA(A0m-3  spectrum  at 
large  m,  in  line  with  observations.  However,  the  situation 
will  generally  be  more  complicated  than  this,  in  which  case 
detailed  knowledge  of  wave  phase  speeds  and  the  meteoro¬ 
logical  environment  is  needed  to  assess  the  wind-shifting 
model  properly.  Such  efforts  are  beyond  the  scope  of  this 
initial  assessment  of  the  process. 

Discussion 

Comparison  With  Existing  Theories 

Saturation  models  cite  attenuated  spectra  as  evidence  that 
the  wave  field  is  unsaturated,  yet  often  there  are  few  clues  as 
to  why  this  should  be  so.  For  example,  unsaturated  spectra 
in  the  stratosphere  are  difficult  to  reconcile  when  the  spectra 
observed  lower  down  in  the  troposphere  are  saturated  [e.g., 
Tsuda  et  al.,  1991]. 

The  wind-shifting  model  supports  the  postulate  that  the 
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wave  field  is  unsaturated  in  such  cases  and  provides  a 
dynamical  explanation  as  to  why  and  where  this  occurs. 

Refractive  downshifting  (j3  =  —  1)  induces  spectral  atten¬ 
uation  by  reducing  wavenumbers  and,  more  importantly  in 
this  case,  attenuating  wave-field  variances.  It  produces 
strong  attenuation  of  the  normalized  shear  variance  a2,  for 
which  the  attenuating  effect  is  approximately  triple  that  for 
u'*  (compare  (4)  and  (5)).  It  is  a2  which  determines  whether 
the  wave  field  (either  one  wave  or  a  collection  of  waves)  is 
saturated  [Desaubies  and  Smith,  1982;  Hines,  1991a],  and  so 
the  reduction  of  a2  by  downshifting  produces  an  unsaturated 
wave  field. 

This  explanation  can  also  be  cast  within  the  Doppler- 
spread  theory  of  Hines  [1991b],  whereby  reductions  in  u'1 
over  the  whole  wave  spectrum  due  to  downshifting  produce 
reduced  spreading  of  wavenumbers  to  dissipative  vertical 
scales  and  reduced  a2  variances,  which  imply  increased 
wave  field  stability. 

Hines  [1993a]  formulated  a  theory,  based  on  Doppler 
spreading,  to  explain  attenuated  spectral  densities  in  the 
stratosphere.  The  idea  depends  upon  the  separation  of  peak 
vertical  wavenumbers  in  the  source  spectrum  (m0)  and  in 
the  Doppler-spread  spectrum  ( mc ,  equivalent  to  m#  in  (8)). 
If  mc  «  m0,  then  the  Doppler-spread  spectra  tended  to  the 
aA(N)m~ 3  shape  at  m  »  mc  for  a  variety  of  source 
spectra  that  he  considered.  However,  when  mc  »  m0  for 
these  same  source  spectra,  the  resulting  Doppler-spread 
spectra  at  large  m  grew  with  altitude,  qualitatively  in  line 
with  stratospheric  observations  [e.g.,  Wilson  et  al.,  1990]. 
Hines  [1993a]  coined  the  latter  effect  “pseudosaturation.” 
Similar  effects  can  occur  within  diffusive  models  of  wave 
spectra  [Weinstock,  1990]. 

Hines  [1993a]  did  not  speculate  on  what  processes  might 
produce  the  necessary  conditions  (i.e.,  mc  5S>  m0)  for 
pseudosaturated  stratospheric  spectra  [e.g.,  Wilson  et  al., 
1991b;  Senft  et  al.,  1993].  The  wind-shifting  model  can 
self-consistently  trigger  or  eliminate  pseudosaturation  con¬ 
ditions  in  line  with  observations. 

When  downshifting  occurs,  wavenumbers  in  the  source 
spectrum  all  decrease,  so  that  the  peak  wavenumber  m0  of 
the  source  spectrum  must  also  decrease.  Wind-shifting  the¬ 
ory  also  predicts  a  reduction  in  the  variance  beneath  the 
source  spectrum,  due  to  EP-flux  conservation.  This  reduc¬ 
tion  in  m7^  leads  to  reduced  Doppler  spreading  and  a 
consequent  increase  in  mc  [Hines,  1991b],  These  oppositely 
directed  wavenumber  shifts  may  then  lead  to  situations 
where  the  environment  changes  from  mc  «  m0  and  unat¬ 
tenuated  spectra  at  large  m  (e.g.,  in  the  troposphere),  to 
mc  5*>  m0  and  attenuated  spectra  at  large  m  (e.g.,  due  to 
strong  downshifting  in  the  stratosphere),  and  thereafter  to 
increases  in  spectral  intensity  with  height  until  limiting 
spectral  shapes  return  when  mc  «  m0  is  reattained  (e.g.,  in 
the  mesosphere).  This  may  also  impact  upon  the  measured 
spectral  index  q,  as  discussed  in  the  next  section. 

Hostetler  and  Gardner  [1994]  measured  attenuated  strato¬ 
spheric  spectra  and  also  noted  an  excess  of  energy  at  small 
wavenumbers  m .  They  argued  that  this  may  reflect  a  gravity 
wave  source  spectrum  dominated  by  small-m  waves  pro¬ 
duced  by  geostrophic  adjustment  of  the  jet  stream,  thus 
accounting  for  their  attenuated  spectra  via  the  source- 
spectrum  arguments  of  Weinstock  [1990]  and  Hines  [1993a], 
The  wind-shifting  model  can  account  for  these  observations 
without  invoking  a  predominant  small-/n  gravity  wave 


source.  Concentration  of  energy  at  small  m  arises  due  to 
refraction  of  vertical  wavenumbers  to  small  values  by  the 
background  stratospheric  winds  (e.g.,  Figure  10),  while 
spectral  attenuation  follows  not  only  from  mc  »  mQ  effects 
but  also  from  reductions  in  variance  owing  to  conservation 
of  EP  flux  under  such  wind-induced  refraction.  Indeed, 
Eckermann  et  al.  [1995]  presented  direct  experimental  evi¬ 
dence  of  refractive  downshifting  of  m  values  in  the  upper 
stratosphere.  Furthermore,  while  the  wind-shifting  model 
has  been  able  to  account  for  some  observed  vertical  and 
seasonal  changes  in  wave  spectra  through  observed  changes 
in  the  background  wind  profile,  source  arguments  must 
invoke  a  time-varying  (and  possibly  even  a  height  varying) 
source  spectrum  to  explain  these  measurements. 

Kuo  and  Lue  [1994]  proposed  a  theory  which  depends 
upon  the  term  w'  dUldz  in  the  linearized  zonal  momentum 
equation.  This  term  is  often  discarded  in  various  approxima¬ 
tions,  including  the  ones  used  here.  On  integrating  general 
forms  of  the  linearized  fluid  equations  appropriate  to  gravity 
waves  in  a  numerical  model,  Huang  et  al.  [1992]  and  Kuo  et 
al.  [1992]  simulated  spectra  which  were  similar  to  observa¬ 
tions.  To  shed  light  on  their  findings  in  regions  of  wind  shear, 
Kuo  et  al.  [1992]  derived  an  analytical  wave  shear  model 
based  on  the  w'  dUldz  interaction  effect  for  a  single  wave, 
in  which  the  shear  could  change  h' 5  with  height  for  a 
constant  w'1  value.  They  found  that  their  model  could 
account  for  their  observed  variations  of  h77  with  height.  Kuo 
and  Lue  [1994]  later  applied  the  idea  to  explain  attenuated 
spectra  at  large  m  observed  with  the  MU  radar  in  the  lower 
stratosphere. 

In  deriving  their  analytical  equations,  Kuo  et  al.  [1992] 
omitted  the  pressure  gradient  term  in  the  zonal  momentum 
equation,  which  is  important.  When  pressure  gradient  terms 
are  retained  in  all  the  momentum  equations,  one  obtains  a 
polarization  relation  between  the  peak  velocity  amplitudes 
u'  and  w'  of  the  form  [e.g.,  Gossard  and  Hooke,  1975] 


u  = 


irn-T  +  ( ch  -  U)  U 
1  -  (c*  -  U)2ICl 
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where  Cs  is  the  speed  of  sound  and  T  =  - 1/2 Hp  +  g/Cj  is 
Eckart’s  coefficient.  The  dispersion  relation,  for  d2/dzz  U  = 
0  and  ignoring  the  Earth’s  rotation,  becomes 
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(12) 


Note  that  retaining  only  the  first  term  in  (12)  yields  (1). 

Equations  (11)  and  (12)  reveal  that  the  effect  of  the  shear 
in  modifying  u'  in  the  presence  of  constant  tv'  is  only 
significant  when  c“‘  =  ch  -  U  is  large  and  m  is  conse¬ 
quently  very  small.  In  fact,  the  effects  of  the  shear  are  very 
small  for  vertical  wavelengths  2  n/m  s  50  km  for  typical 
values  of  background  atmospheric  variables  and  using  the 
largest  dUldz  values  quoted  by  Kuo  et  al.  [1992].  In  light  of 
this  analysis,  there  seems  to  be  no  direct  way  that  shearing 
effects  such  as  these  can  significantly  influence  the  intensity 
of  the  spectrum  at  the  larger  vertical  wavenumbers,  as 
argued  by  Kuo  and  Lue  [1994]. 

Spectral  Indices  of  Attenuated  Spectra 

A  major  weakness  of  the  current  spectral  wind-shifting 
theory  is  that  while  it  predicts  the  occurrence  of  attenuated 
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spectra,  it  has  proved  less  successful  in  estimating  then- 
spectral  indices  q.  The  model  of  Fritts  and  Lu  [1993] 
predicts  increases  in  q  to  values  exceeding  3  (see  Figure  3). 
While  this  is  evident  in  some  spectral  data  [e.g.,  Tsuda  et  al., 
1994b],  more  usually  nominal  q  values  —2-2.5  are  observed 
[e.g.,  Tsuda  et  al.,  1991;  Marsh  et  al.,  1991;  Senft  et  al., 
1993;  Hostetler  and  Gardner,  1994;  Mitchell  et  al.,  1994]. 

Fritts  and  Lu  [1993]  developed  their  model  with  a  view  to 
parameterizing  gravity  wave  drag  within  general  circulation 
models  and  thus  did  not  formulate  it  with  the  analysis  of 
observational  spectra  in  mind.  Wind-shifting  effects  were 
incorporated  in  a  simple  way,  whereby  the  spectral  density 
at  each  harmonic  m,  was  associated  with  a  wave  or  waves  of 
wavenumber  mh  which  may  not  be  appropriate.  For  exam¬ 
ple,  in  the  Doppler-spreading  model,  the  spectrum  at  large  m 
is  formed  by  statistical  effects  and  cannot  be  associated  with 
discrete  waves  at  each  /n,  [Allen  and  Joseph,  1989;  Hines, 
1991b], 

However,  while  further  observations  are  needed  to  char¬ 
acterize  these  q  data  better,  the  data  to  date  certainly 
suggest  strongly  that  in  addition  to  the  attenuation  of  spec¬ 
tral  density  at  large  m,  downshifting  frequently  yields  a 
differently  shaped  spectrum  as  well  (i.e.,  q  ¥=  3).  We  shall 
discuss  some  possibilities  as  to  how  and  why  shape  changes 
may  occur,  particularly  the  occurrence  of  q  <  3  which  the 
wind-shifting  model  employed  here  does  not  predict. 

As  argued  earlier,  reductions  in  u' 1  due  to  downshifting 
will  increase  the  cutoff  wavenumber  mc  according  to  the 
Doppler-spreading  model  (or,  equivalently,  increase  m* 
according  to  the  saturation  model)  [Smith  et  al.,  1987;  Hines, 
1991b].  This  wavenumber  (mc  or  /w*)  marks  a  transition 
zone  where  the  spectral  shape  changes,  from  linear  variabil¬ 
ity  with  matm«  mc  or  m*,  toward  the  aA{N)m  ~3  shape 
at  m  :»  mc,  or  m*  [Hines,  1991b;  see  also  Figure  3].  Thus 
an  increase  inmc  due  to  the  attenuating  effects  of  down¬ 
shifting  acts  to  move  this  spectral  transition  zone  to  larger 
wavenumbers.  This  may  explain  observations  of  q  <  3  at 
large  m  for  attenuated  spectra,  since  they  may  now  lie  within 
the  transition  zone  where  - 1  <  q  <  3.  This  idea  differs  from 
the  Fritts  and  Lu  [1993]  model  where  downshifting  also 
“downshifts”  m*,  which  is  appropriate  for  instantaneous 
wind  shifts  but  may  not  persist  once  waves  grow  in  ampli¬ 
tude  and  interact  with  each  other  [e.g.,  Hines,  1991b]. 

Extratropical  gravity  wave  data  from  the  upper  strato¬ 
sphere  provide  some  support  for  this  general  scenario. 
Observations  reveal  an  annual  cycle  in  wave  variances,  with 
largest  variances  in  winter  and  smallest  values  in  summer 
[Hirota,  1984;  Mitchell  et  al.,  1991;  Wilson  et  al.,  1991b; 
Eckermann  et  al.,  1995],  As  downshifting  occurs  in  both 
summer  and  winter  (see  Figure  10),  wind  shifting  does  not 
cause  these  seasonal  changes  in  the  variance.  Instead,  they 
occur  mainly  through  seasonal  variations  in  the  exp  (J 
dzlHp)  term  [Eckermann,  1995].  This  annual  variation  of 
wave  variance  should  produce  smaller  mc  (m*)  values  in 
winter  than  in  summer,  and  thus  the  large-m  portion  of  the 
spectrum  should  be  closest  to  the  limiting  m  “3  shape  in 
winter,  when  the  inequality  m  55>  mc  is  most  likely  to  be 
satisfied.  Consequently,  q  should  be  closest  to  (farthest 
from)  3  in  winter  (summer).  Upper  stratospheric  observa¬ 
tions  do  indeed  reveal  an  annual  cycle  in  q,  with  q  nearest  3 
in  winter  and  q  —  2  during  summer  [Eckermann,  1990b; 
Marsh  et  al.,  1991].  Given  that  q  — »  3  as  variances  increase 
into  the  mesosphere,  as  discussed  earlier,  then  there  appears 


to  be  evidence  of  a  correlation  between  reduced  (increased) 
variances  and  q  values  smaller  than  (close  to)  3,  which  is 
somewhat  consistent  with  associated  shifts  in  mc  (m*). 

Thus  even  though  the  seasonal  variance  asymmetry  dis¬ 
cussed  above  is  not  caused  by  wind  shifting,  attenuation  of 
variance  by  wind  shifting  can  produce  similar  sorts  of  effects 
and  thus  may  account  for  the  q  <  3  data.  These  arguments 
assume  that  a  characteristic  m  _3  spectrum  (either  saturated 
or  pseudosaturated)  always  returns  at  m  »  mc,  in  line  with 
current  theories  [e.g.,  Smith  et  al.,  1987;  Allen  and  Joseph, 
1989;  Hines,  1991b,  1993a],  Yet  some  of  the  observations 
suggest  statistically  significant  measurements  of  q  ¥=  3  at  m 
values  much  larger  than  predicted  mc  values  [e.g.,  Wu  and 
Widdel,  1990;  Tsuda  et  al.,  1991], 

We  have  only  considered  waves  propagating  parallel  to 
the  background  wind,  whereas  waves  will  exist  at  various 
angles  <p  with  respect  to  the  background  flow  (see  Figure  la). 
More  complete  theoretical  studies  of  three-dimensional 
propagation  of  a  wave  spectrum  within  different  types  of 
background  wind  profiles  will  be  needed  to  anticipate  the  full 
effects  of  these  wind-shifting  processes  on  spectral  densities 
and  shapes  [e.g.,  Muller  et  al.,  1986;  Huang  et  al.,  1992; 
Sato  and  Yamada,  1994].  The  wind-induced  changes  in 
variance  with  height  (e.g..  Figure  7)  imply  a  highly  nonsta¬ 
tionary  data  series,  which  can  also  produce  important 
changes  in  spectral  shapes  [Eckermann,  1990a],  as  can  other 
data-dependent  effects  such  as  time  averaging  and  aliasing 
[e.g.,  Lefrere  and  Sidi,  1990;  Gardner  et  al.,  1993]. 

Summary  and  Conclusions 

Here  it  has  been  proposed  that  the  occurrence  of  attenu¬ 
ated  vertical-wavenumber  spectral  densities  at  large  m  in  the 
observational  literature  is  caused  by  the  specific  nature  of 
the  vertical  profile  of  background  wind  velocity  within  and 
beneath  the  observed  region.  When  attenuated  spectra  oc¬ 
cur,  the  background  winds  U  are  such  that  certain  waves  are 
removed  from  the  spectrum  through  critical-level  effects 
lower  in  the  atmosphere.  This  leaves  a  directional  wave  field 
higher  in  the  atmosphere  for  which  either  ch  <  U  or  ch  >  U 
is  preferred. 

Attenuated  spectra  are  observed  when  the  shear  is  such 
that  U  “moves  away”  from  ch  as  the  waves  propagate 
upward,  and  so  vertical  wavenumbers  m  ~  N/\ch  —  U  cos 
<p\  decrease.  This  is  the  “)3  =  —1”  (downshifting)  scenario, 
where  /3  is  defined  by  (7).  On  the  other  hand,  little  spectral 
attenuation  and  shape  variability  is  observed  in  the  j8  =  + 1 
(upshifting)  case,  where  m  values  increase  with  height  due  to 
U  moving  closer  to  phase  speeds  ch  as  the  waves  propagate 
upward. 

Theoretically,  changes  in  U  which  yield  /3  =  —  1  also 
produce  attenuated  wave  field  variances  (to  conserve  EP 
flux)  and  reduced  vertical  wavenumbers  (due  to  refraction). 
Application  of  these  ideas  to  spectral  models  leads  to  reduc¬ 
tions  in  spectral  densities  and  changes  in  spectral  shape  at 
large  vertical  wavenumbers  m  [Fritts  and  Lu,  1993]  which 
are  in  qualitative  accord  with  observations.  When  /3  =  +1, 
variances  and  wavenumbers  increase,  triggering  instabilities 
which  lead  to  dissipation  of  wave  variances  back  to  margin¬ 
ally  stable  (saturated)  values.  In  such  cases,  unattenuated 
spectral  shapes  at  large  m  are  expected  and  tend  to  be  found 
in  observational  data. 

Consequently,  the  wind-shifting  theory  outlined  here, 
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while  admittedly  simplistic,  nevertheless  succeeds  in  mod¬ 
eling  the  occurrence  and  qualitative  form  of  vertical  wave- 
number  spectra  at  large  m  observed  within  different  back¬ 
ground  environments.  These  findings,  if  taken  as  valid,  can 
be  viewed  as  initial  observational  support  for  efforts  to 
incorporate  wind-shifting  effects  into  spectral  parameteriza- 
tions  of  gravity  wave  variability  within  the  middle  atmo¬ 
sphere  [Fritts  and  Lu,  1993].  Indeed,  these  wind-shifting 
processes  are  spectral  analogs  of  the  refraction  and  dissipa¬ 
tion  effects  for  single  waves  considered  by  Lindzen  [1985]. 
However,  the  precise  changes  in  spectral  shape  at  large  m 
are  not  always  well  predicted  by  the  current  models  of 
wind-shifted  spectra  (e.g.,  the  spectral  index  q).  Therefore 
further  observational  and  theoretical  studies  are  needed  to 
fully  gauge  the  extent  to  which  wind-shifting  ideas  can 
explain  observed  spectral  variability  at  large  m. 
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